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<3 ! Abstract 

We theoretically investigate breathing oscillations of fermion condensates in anisotropic har- 
monic oscillator traps. When anisotropy of the trap is not large, the Fermi gases behave as three- 
dimensional (3D) gases, and those breathing oscillations are studied in the previous work l[; when 
the trap is super-anisotropic, the fermions show quasi-low-dimensional (QLD) properties differently 
from the 3D gases because of variation of one-particle energy-level densities, which exhibit specific 
types of shell-structures, QLD shell- structures, in super-anisotropic systems. In the present work, 

> 

oo : we study the effects of the QLD shell-structures on the breathing oscillations of the week-interacting 

fermion condensates with same number densities. The QLD shell-structures are described in the 

"xT ' 

semi-classical Thomas-Fermi approximation (TFA) extended to the super-anisotropic systems, and 

o . 

ON . the collective frequencies of the breathing oscillations are estimated in the scaling and sum-rule 

o 

methods with the perturbation theory and extended TFA up to the first order in a coupling con- 



stant of inter-atomic interaction. As a result, we reveal that the effects of the QLD shell-structures 
on the collective frequencies disappear in transverse modes in both of the oblate and prolate de- 
formed traps, and appear only in longitudinal modes with hierarchic structures associated with the 
trap anisotropy. We also demonstrate time-evolution of the oscillations in the present framework. 
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I. INTRODUCTION 



Development of the trapping and cooling techniques of atoms yields many and various 
studies on trapped ultracold atomic gases. One can now study the Bose- Einstein condensates 



(becs) lalaJlL* 



of them 



0, fl E 



ie degenerate Fermi gases (or fermion condensates) [5|, |6| , and mixtures 



101 ]. In particular, the two-component Fermi gases are recently well- 



studied in various interaction regimes: In attraction regime, the Feshbach resonance method 
realize the atomic pair condensates in the BCS theory and molecular BEC (the BCS-BEC 



crossover 



12|; in strong repnlsion reg nne, the repnlsion indnce a phase transition from 
paramagnetic states to ferromagnetic and phase-separation states M- In addition, the 

D) atomic gases are also realized in super-anisotropic traps [141 . 



quasi- 
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ow- dimensional 



13, 
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221 ]. which often show quite different behaviors from the three- 



dimensional (3D) gases. Those systems offer great infrastructure of fundamental study on 
quantum many-particle systems. 

One of important features in the quantum many-particle systems is the collective excita- 
tions, which are often sensitive to interaction and structure of microscopic quantum states, 
and must reflect the QLD properties. In particular, the breathing oscillations give impor- 
tant diagnostic signals for global properties of the trapped quantum gases in the actual 
experiments. Those collective excitations are theoretically treated with the time-dependent 
mean-field theory, which is reduced to the random-phase approximation (RPA) for small 



amplitude excitations 23j, |24j, |25|]. When the one-particle spectrum is regular and exhibits 



the collective excitations in the long-wavelength range, one can use simpler methods to cal- 

261 ] and the 



culate the frequencies of the collective oscillations: e.g. the sum-rule method 



scaling met 
In Ref. 



rod 



291 ]. one of the authors and his collaborator study spin excitations on dipole 
and monopole oscillations in the 3D two-component Fermi gases in the spherical trap using 
the scaling method. In that work, we discuss the mixing of the in-phase and out-of-phase 
oscillation modes and show the relations between the oscillation frequencies and the phase 



of the two-component Fermi gases [131 ] . 

In Ref. [l], furthermore, we study the breathing oscillations of the two-component fermion 
condensates in anisotropic harmonic oscillator traps with various trap frequencies, and 
demonstrate the mixing behavior of breathing oscillations in the longitudinal and trans- 
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verse directions. In that work, however, we assume that anisotropy of the traps is not so 
large, and apply the conventional Thomas-Fermi approximation (TFA) for the 3D gases to 
predict the collective frequencies of the breathing oscillations in the scaling method. Thus 
it is needed to improve the approach when we study the collective oscillations in the QLD 
systems. 

In this paper, then, we extend those studies on the breathing oscillations of the 3D 
gases to the super-anisotropic deformed gases with the QLD properties, which gradually 
appear as the trap anisotropy increases. The QLD properties should originally be given in 
the microscopic approach, where the quantum many-body properties can be described with 
the shell-structures of the one-particle spectrum by using the one-particle picture, so that, 
in principle, we can consider the QLD properties as exhibition of specific types of shell- 
structures, QLD shell-structures, in the super-anisotropic systems. The precise properties 
of the QLD shell-structures are studied and discussed in the present work. 

The aim of the present work is to reveal the QLD properties on the breathing oscillations 
of the fermion condensates and describe the crossover behaviors between the 3D and QLD 
gases. To simplify the subject, we focus on the breathing oscillations of the week-interacting 
fermions with same number densities, i.e. the in-phase breathing oscillations, and use the 
perturbation theory and semi-classical TFA to calculate the parameters needed in the scaling 
and sum-rule methods. 

The content of this paper is as follows. In Section [Til we gi ye the theoretical framework in 
this work: theoretical explanations of the trapped Fermi gases and breathing oscillations and 
formulation of the scaling and sum-rule methods. In Section lTTTI we consider the ground state 
properties reflecting the QLD shell-structures by developing TFA to the super-anisotropic 
systems in order to calculate the parameters needed in the scaling and sum-rule methods. 
In Section IIVt we show calculational results of the collective frequencies and demonstrate 
time-evolution in the present framework. In Section |V], we give summary and outlook. 
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II. THEORETICAL FRAMEWORK 



A. Trapped Fermi gases 

Let us take two-component gases of fermionic atoms of Ax and A 2 species with same 
atomic masses, mi = m 2 = m, in same trap potentials, V\{r) = V 2 {r) = V(r). The 
Hamiltonian is denoted as 

H W =Y^jdr «rt ( r ) ^_ 1 v 2 + V(r)) ¥ a (r) + H int , (1) 

where (r) is the fermion field operators of atoms, and we choose units for the 
Plank constant h and atomic mass m as h = 1 and m = 1. At ultra-low temperature, the 
inter-atomic interaction part H int in Eq. ([I]) becomes 

= <7 /^ ^(^^(r)^^)*!^) (2) 

with the pseudo-potential of the low-energy s-wave scattering. Here we determine the cou- 
pling constant g in Eq. (j2|) with the s-wave scattering length an as g = Anxayi, and neglect 
the inter-atomic interactions between same kinds of atoms because of the Pauli blocking 
effect. The trap potential V(r) in Eq. ([T]) is a cylindrical harmonic oscillator potential, 

V(r) = ^y c +uy z ), (3) 

where we introduce the Cartesian coordinates, r = (r x , r y , r z ), and r 2 c = + r^, and assume 
the axial symmetry around the z-axis. In this paper, we call longitudinal and transverse 
directions for the parallel and orthogonal directions to the z-axis, respectively. Ratio of the 
trap frequencies, uj c /u z , decides anisotropy of this system. 

To simplify the model, we consider the symmetric two-component gases with same atomic 
densities, pi(r, t) = pi (r, t) = p(r,t), and use the time-dependent Hartree-Fock (HF) mean- 
field approximation [30j . In the HF approximation, the two-body interaction in Eq. can 
be rewritten into a one-body interaction with the self-consistent mean-fields of the atomic 
densities, and the many-body effects are included in the mean-fields as mean-values. Then 
we can omit the subscripts for the atomic species by introducing a redefined system for the 
each component gas: 

H = J dr *t( r ) ^_I V 2 + V(r)^j (r) + H int (4) 



with 

H int = g J dr U\r)p{r,t)^{r) - ^(p(r,t)) 2 ) , (5) 

where the atomic density, p(r,t) = (ty^ty}, must be determined self-consistently. Here we 
should comment that symmetric /-component gases, / is an integer, take the same formula- 
tions in this paper with a renormalized coupling constant, g —> (I — l)g; thus the results in 
this paper can apply to any symmetric /-component gases. 

In the HF approximation (Eqs. fll]) and (jSJ), the ground state |<3> ) is given by the Slater 
determinant 

l $ o) = (n4)kac> (6) 

n<?ip 

with the fermion creation operator 

C t = J dr ^(r)0„(r), (7) 
where the normalized one-particle wave functions 4> n {r) are determined by the HF equations, 

-~ V 2 + V(r) + gp(r)) n (r) = e n n (r), (8) 
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where p(r) is the atomic density in the ground state. The Fermi level rip in Eq. ([6]) decide 
the atom number 

N = j drp(r)= I drJ2\Mr)\ 2 . (9) 

n<np 

According to Eq. (jSJ), the ground energy E s = ($o\H |^o) becomes 

E S = [ dr J2 <f>* n (r) (-^V 2 + V(r)) Mr) + E int (10) 

J n<n F ^ ' 

with the interaction energy 

E int = (<5> \H int \%) = \jdr (p(r)) 2 . (11) 
The HF approximation is valid when the gases are weekly interacting and dilute, g ~ 

— 1/2 

and |a 12 | <C uj c ,z , as assumed in the present work. If the gases have strong correlations, 
e.g. the pair correlation in the BCS-BEC crossover Q ^ the QLD colons in the 
non-weekly interacting QLD systems 



221 ] . one should use another different method including 
the correlation effects. Furthermore, if the gases go to dense by deformation, the pseudo- 
potential in Eq. (jSJ) for the 3D scattering may not be available because of the QLD scattering 



effect. In this case, one should consider more realistic inter-atomic interactions instead of 
the pseudo-potential in Eq. (j2J). Those detailed discussions are beyond the purpose of this 
work, and omitted in this paper. 

B. Breathing oscillations 

The breathing oscillations are originally defined as compressive oscillations of fluids; espe- 
cially, in spherical systems, they mean the monopole oscillations in the sense of the multipole 
expansion. However, in anisotropic systems, the monopole oscillations are inseparably mixed 
with incompressive oscillations (e.g. the quadrupole oscillations) an d we need to redefine 
the breathing oscillations as the mixtures of those oscillations as taken in this paper. 

In the present work, we assume cylindrical deformation in excited gases according to 
the cylindrical trap potential in Eq. (J3]). Then the breathing oscillations appear in the 
expectation values of projected mean square radius operators 

R c>z = J dr &(r)f*Mr)- (12) 

In the actual experiments, those expectation values are directly observed with the absorption 
imaging method, and the excited gases are artificially generated with sudden change of the 
trap potentials. 

When the variations of the expectation values of the R c>z in Eq. (fT2|) , 

s c>z = ($(t)\R c>z \$(t)) - ($ \R C>Z \%), (13) 

exhibit small amplitude oscillations, i.e. when the system shows the minimal breathing 
oscillations, we can describe these oscillations as 

Sc,z = J2 C »(®v\ R cA®») e ~ iEvt + C - C - ( 14 ) 
v 

with the time-evolving excited state 

|$(t)) = |$o) + ^a|^)e-^, (15) 

V 

where we introduce the excited states \Q U ) with the excited energy E u and the minimal 
amplitudes C v determined by the initial condition. As shown in Eq. (Tl^I) . those oscillations 
directly reflect the information of the ground and excited states. The oscillations of the s c (t) 
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and s z (t) in Eq. tHM are called transverse and longitudinal oscillations, respectively, which 
correspond to the monopole oscillations when s c (t) = s z (t) and the quadrupole oscillations 
when s c (t) = —s z (t)/2. In general, the transverse and longitudinal oscillations are coupled, 
so that one should consider the normal mode oscillations given by the normal mode operators 
R c , z , which are described as linear combinations of R c z and asymptotically correspond to 
i? C)2 at a decoupled limit. In the collective oscillations, the transition strength functions 
($o \Rc,z p&v) are localized in a small range of the excited energy; then the corresponding 
normal mode oscillations behave like harmonic oscillations according to Eq. (CHI) , and one 
can theoretically estimate the frequencies by using the scaling and sum-rule methods. 

C. Scaling and sum-rule methods 

We now formulate the scaling and sum-rule methods for the breathing oscillations of the 
fermion condensates, which are defined as the mixed oscillations of the s C)2 in Eq. ( fl3l) as 
explained above. 

According to Eq. (Tl3|) . the behaviors of the breathing oscillations are determined from 
the time-dependent excited state | In the scaling method, the excited state is given 

by the scale transformation of the ground state |$o) : 

|$(t)) = e -^C e -MH,R c ] + X z [H,R z ])/2^ (16) 

with the projected mean square radius operators R CyZ defined in Eq. f fl2|) and the Galilei 
transformation factor 

where A Cj2 represents the time derivative of A Cj2 . The scale parameters \ c ,z(t) in Eq. (fl6|) 
indicate the collective coordinates of the transverse and longitudinal oscillations, and are 
proportional to the projected mean square radius s c>z in Eq. (fl3l) as a result. 

In the HF approximation, the ground state is given by the Slater determinant in Eq. ([HD; 
then the excited state in Eq. ( fT6l) becomes 

i^))=(n^))i vac > 

ra<jip 

with 

4(t) = j drV\r)4> n (r,t), (19) 



where we introduce the scaled wave functions defined as 

4> n (r, t) EE e -^+A^)/2 e A c +A,/ 2( ^ ^A Cr , c) ^ J (2Q) 

with r c ee (r x ,r y ). These scaled wave functions satisfy the normalization and continuity 
conditions. 

The time- dependent variational principle, 

5 J dt{${t)\ij t - H\${t)) =S J dt £[A,A] =0, (21) 

gives the equations of motion of the scale parameters \ c ,z(f) in Eq. ffTB]) . According to 
Eq. (118jl . the variational space in Eq. (1211) is limited in the configuration space of the Slater 
determinants, so that the scaling method is in the time-dependent HF theory. 

By substituting Eq. ([TBI into Eq. (l2Tj) . the Lagrangian in Eq. (J2TJ) can be written as 

£[A, A] = - H, (22) 



J=C,Z 



where we define the mass parameters 



with 



M W = e -^i M . (23) 



Mj ee J drh 2 jP {r) (24) 

and the excited energy 

H[\X\ = (m\H\m)-E g (25) 
with the ground energy E g in Eq. ( flOl) . The excited energy in Eq. ( |25l) is calculated as 

(m\H\Ht)) = lY, M J X) ^ +v ^ ( 26 ) 
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with the potential parameter 

V[A] = ( e2XjR J + e_2Aj ^') + e2Ac+Az ^nt, (27) 

where Kj and £/j are the kinetic and trap potential energies of the ground state in the 
j-direction motion, 

Kj ee Jdr ($ |* f (r) *(r)|$ ) (28) 
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and 

Uj = Jdr <$o|^ f (r) Q^l) *W| $ o>- (29) 

In the present work, we study the minimal oscillations (X c ,z C 1). To treat the minimal 
oscillations, we expand the excited energy in Eq. (1251) by A Cj2 and take the expanded terms 
up to the second order: 

H[\, A] ~ \>?B\ + l* T C\ (30) 

with A = (A c , A^) T , where the superscript T indicates the transposition. The B and C in 
Eq. (|30|) indicate the mass parameters, 



(31) 




and the restoring forces, 

c ^ f 4(K c + U c )+AE int 2E int 

\ 2E- mt ±{K Z + U z ) + E int , 

respectively. Here one should note that the mass parameters B in Eq. (13T1) is diagonal 
because [ [H, i? c ] , R z \ = and [ [H, R z \ , i? c ] = 0; in addition, the first order terms of A in 
Eq. (130]) must vanish owing to the virial theorem, which derives 

K C ~Y1 J dr P( r )f^ r ) = K C -U C + E int = (33) 



and 



with the virials 



K z - Jdr p(r)f z (r) = K z - U z + ^E int = (34) 

fi(r) = ^^r(yir)+gp(r)). (35) 

In order to diagonalize the 7i[A, A] in Eq. ( l30i) . we introduce a complete set of orthonormal 
eigenvectors, 

{B- 1 ' 2 CB- 1 ' 2 )r 1j = nfa (36) 
for j = c and z with the eigenvalues Q CtZ . Then the A in Eq. (|30|) can be described as 

A = B-^ Yl Wi ( 37 ) 
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with 

Uj = t]jB 1/2 \. (38) 

By substituting Eq. fl37|) into Eq. (130]) . the TC in Eq. fl30l) can be written as the diagonal 
form, 

ft = Ia t A + ^A T (7A (39) 
with the normal mode vector A = (u c ,u z ) T and the eigenvalue matrix 

~ (n 2 c \ 

C=\ . (40) 

According to Eq. fl39l) . the f2 C)2; in Eq. (136]) correspond to the collective frequencies of the 
normal modes. Using the normal mode representation, we obtain the decoupled scale trans- 
formation 

= e -i(^cRc+u z R z )/2 e -(u c [H,R c ]+u z [H,R z ])/2^^ ^ 

with the normal mode operators 

R 5 = r)jB- 1/2 R, (42) 

where R = {R c , R Z ) T . 

In principle, the scaling and sum-rule methods predict the same frequencies of the oscil- 
lations as shown in Appendix |A] However the scaling method can describe an approximated 
time-evolution of the oscillations as shown in Eq. (1TB]) . and also give coupled oscillations 
with several modes, i.e. determination of the Rj in Eq. f)42p . differently from the sum- rule 
method; on the other hand, the sum-rule method can describe more detailed information 
of the oscillations, e.g. collectivity of the oscillations, with additional data of the moments, 
differently from the scaling method. Thus the scaling and sum-rule methods have comple- 
mentary roles to study the collective oscillations. 

According to the above formulation of the scaling and sum-rule methods, behaviors of the 
breathing oscillations depend only on the parameters, A4. c ,z, K C Z) U C)Z , and E- m t, determined 
from the ground state properties. In fact, when these parameters are given, one can calculate 
the collective frequencies in Eq. (|36l) and predict the precise behaviors via the time- dependent 
excited state in Eq. f|T6|) by solving the dynamical system in Eqs. ff22|) and fl30|) . 
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III. GROUND STATE PROPERTIES 

Let us calculate the parameters needed to predict the behaviors of the breathing os- 
cillations in the scaling and sum-rule methods. These parameters are determined only 
from the ground state properties as explained in the previous section, and reflect the QLD 
shell-structures through the ground state properties. Thus we here consider the ground 
state properties and their QLD shell-structures by developing the semi-classical TFA to the 
super-anisotropic systems, and then determine the parameters by using this approach. 

A. QLD shell-structures and extended TFA 

The QLD shell-structures are defined as specific types of the shell-structures in the super- 
anisotropic systems. In general, the shell-structures can be well-defined in some asymptotic 
states, and it is useful for theoretical study to introduce such asymptotic states. Thus we 
here introduce the asymptotic states to study the QLD shell-structures. In the present 
work, we choose the non-interacting states (g = 0) for the asymptotic states and use the 
perturbation theory to include the inter-atomic interaction. The non-interacting asymptotic 
states implicitly have asymptotic variable-separation and corresponding asymptotic quan- 
tum numbers, so that the QLD shell-structures must be well-defined in the non-interacting 
asymptotic states and produce the QLD properties at the super-anisotropic limits in this 
framework. 

According to the first order perturbation theory, all physical values are evaluated from 
the non-perturbative state, i.e. the ground state of the non-interacting gases (g = 0). In 
the non-interacting gases, the one-particle energies in Eq. (jSJ) become 



with n = {n x ,n y ,n z } and the corresponding one-particle wave-functions 

n (r) take the variable separation form: 




(43) 




(44) 



j=x,y,z 




(45) 
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with the Hermite functions H n (x). Here it should be noted that the asymptotic variable- 
separation and asymptotic quantum numbers are introduced in Eqs. ( 144|) and fj4"3]) . respec- 
tively. 

In order to develop TFA, we here introduce the Wigner function defined as 

f(r, p)^Jds e-*P s (*t ( r _£)*( P + |)). (46) 
Then the expectation value of a one-body operator 

A = J dr ^ (r)a(f,p)\& (r) (47) 

becomes 

(A) = j dr l dp ^ d) (r lP )f(r lP ) (48) 

with the Weyl classical function a^ cl \r,p) corresponding to the quantum mechanical oper- 
ator a(r,p). In the non-interacting gases, the Wigner function becomes 

f(r,p)=J2He n ) J] fA'-rl'r'j) (49) 

n j=x,y,z 

according to Eqs. (l4*3l) and (jSJ), where we define 

fA'- rl'r tj) = J ds e ' iPjS( Pj ( r i - |) ^ ( r J + ( 5 °) 

and use the zero-temperature Fermi distribution function JF(e n ) = 6 (eF — e n ) in Eq. (149]) 
with the Fermi energy ep and the Heaviside's step function G (x). The Wigner function in 
Eq. (1491) is simplified in TFA by smoothing the energy-level density. 



1. Conventional TFA 

The conventional 3DTFA, explained in Appendix [Bj is valid only when the anisotropy 
of the trap is not so large; otherwise 3DTFA loses the validity owing to the QLD shell- 
structures. In fact, this proposition is clearly demonstrated in the one-particle energy- level 
densities 3l| . We now show them in the oblate and prolate deformed traps in order. 

First we consider the oblate deformed potentials (u z > u c ). Here we take the anisotropy 
a z = u) z /oj c to be an integer to simplify expressions. At that time, the energy-level density 
can be defined as the mean value of the number of energy-levels in an energy width uo c 
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associated with an energy eigenvalue e q = u c q c + eo, where we introduce the quantum 
number q c for the c-direction motion and the zero point energy eo = uj c + u> z /2. It becomes 

n / \ + l)(2g e + 2 - a z q z ) 

= 2^ (51) 

with the maximum quantum number q z = int((e 9 — eo)/u z ) for the z-direction motion, where 
the round-down function int(x) returns the integer part of x. 

By smoothing the energy level density in Eq. ( 15T1) with q c<z pa (e — eo)/uj CjZ , we obtain 

D(s)«^ 3 )(8)^-^-e( £ ), (52) 

ZOu UJ Z 

where -D(e) agrees with that in 3DTFA (see Eq. (1B7j) ). According to the smoothing method, 
this agreement is approved only when the energy level density in e q ^> uj z > u c principally 
contribute to the system, i.e. the anisotropy of the trap is not large; otherwise the QLD 
shell-structures have non- negligible influence. 

In fact, when we consider the super-anisotropic limit, the energy level density at q z = 
(uj z > e q ^> uj c ) principally contribute to the system; then, by smoothing the energy level 
density in Eq. fl5T|) only with q c pa (e — eo)/cu c , we obtain 

(53) 

where D(e) agrees with that in 2DTFA (see Eq. (jB~7|) ). 

In Fig. [U we show the energy level densities when oj z /uj c = 21. When e ^> u z , D(e q ) (the 
open circles) approximately agrees with D^(e) (the dashed line) as described in Eq. (|52[) ; 
otherwise they have some difference reflecting the QLD shell-structures. When e < 31.5 u c = 
(3/2)u z (q z = 0), D(e q ) (the open circles) agrees with D^ 2 \e) (the solid line) as described 
in Eq. (|53|). 

Second we consider the prolate deformed potentials (u c > uj z ). Here we also take the 
anisotropy a c = oj c /uo z to be an integer. At that time, we can use the same way in the oblate 
deformed potentials with exchange of the subscripts (c z) to obtain the energy-level 
density; it becomes 

DM . fe + 1 2 'fe + 2 > . (54) 

According to the same way in the oblate deformed potentials, the energy level density 
in Eq. agrees with that in 3DTFA, D^(e), when the anisotropy of the trap is not so 
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FIG. 1: (color online) The energy level densities when lo z /lo c = 21. The open circles represent 
the values of D(e q ); the dashed and solid lines correspond to the results in 3DTFA and 2DTFA, 
respectively. 

large; otherwise they disagree owing to the QLD shell-structures. Especially, at the super- 
anisotropic limit (q c = and uj c > e q u> z ), the energy level density in Eq. (1541) can be 
written as 

D(e) « D {1) (e) = — 6 (e - cu c ) (55) 

by smoothing q z « (e — eo)/uj z , where -D(e) agrees with that in 1DTFA (see Eq. (1B7I) ). 

In Fig. [21 we show the energy level density when uj c /uj z = 21. When e ^> u; c , D(e q ) 
(the open circles) approximately agrees with D^(e) (the dashed line) as same as that in 
the oblate deformed potentials. When e < 42 u z = 2u c (q c = 0), D(e q ) (the open circles) 
agrees with D^'(e) (the solid line) as described in Eq. fl55l) . The appearance of the QLD 
shell-structures in the prolate deformed systems is clearer than that in the oblate deformed 
systems. 

Here it should be noted that the QLD shell-structures exhibit the hierarchic structures 
associated with the maximum quantum numbers for the narrow direction motions, q z and 
q c , in the both cases of the oblate and prolate deformed potentials as shown in Figs. [1] and 
El respectively. The appearance of those hierarchic structures is a specific feature of the 
QLD shell-structures, and also shown in the following results. In general, we can define 
the QLD shell-structures as the hierarchic structures in one-particle energy-level densities 
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FIG. 2: (color online) The energy level densities when lu c /uj z = 21. The open circles represent 
the values of D(e q ); the dashed and solid lines correspond to the results in 3DTFA and 1DTFA, 
respectively. 

in asymptotic states, where one-particle wave-functions can be denoted as product of two 
wave-functions with different energy-levels and asymptotic quantum numbers (i.e. fine and 



coarse modes) [35]. 



2. Extended TFA 



As shown in the energy-level densities, the conventional 3DTFA should not be applied to 
the super-anisotropic systems owing to the QLD shell-structures. We here introduce a new 
approach to develop TFA to the super-anisotropic systems. 

In order to take into account the QLD shell-structures, we smooth the Wigner function 
in Eq. (H91) only for the c-direction in the oblate deformed potentials and for the ^-direction 
in the prolate deformed potentials: 

f(r,p)= / deGW(e) (56) 
Jo 

with 

G W( e ) = T (e + eSVcpj) ^5 (e - e z ) f z (r z ,p z ; e z ) (57) 
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in the oblate deformed potentials (j = 2) and 

(&\e) = r(e + e%\r,,p,j) £ 5 (e - e x - e y ) \[ ./}(/> /'r < ,) (58) 

n x ,n y j=x,y 

in the prolate deformed potentials (j = 1), where we introduce the two and one dimensional 
classical energies, e^\r c ,p c ) = p\j2 + cu^r^/2 and (r z ,p z ) = p 2 z /2 + ou 2 r 2 /2. 

Here it should be noted that the Wigner function in Eq. (1561) agrees with that in 3DTFA, 

/( 3 )(r,p)^(e cl (r,p)) (59) 

with e c i(r,p) = p 2 /2 + V(r), when e ^> u> c , u z ; thus we can rewrite the Wigner functions in 
Eq. (J5SD as 

f(r,p) = f C de G®(e) + / (3) (r,p)0 (e cl (r,p) - e c ) (60) 
Jo 

in the oblate (j = 2) and prolate (j = 1) deformed potentials with a large cutoff energy e c . 
In this paper, we call the rewritten form of the Wigner function in Eq. (1601) as extended TFA 
(ETFA), which agrees with the original Wigner function in Eq. (1561) when e c > e?. 

The hierarchic structures can also be seen in Eq. (1601 through Eqs. (|57|) and (158]) . so that 
ETFA must reproduce the QLD shell-structures below the cutoff energy; above the cutoff 
energy, ETFA have no shell-structure as shown in Eq. (1591) . However, pragmatically, it is 
no need to treat the QLD shell-structures in all range of the energy-level density because 
the QLD shell-structures have non- negligible influence only in the low-energy range. In fact, 
the results in ETFA with finite cutoff energies show good agreement with those in the exact 
calculations (e c — > oo) as demonstrated in the next section. 



3. Fermi energy and, crossover behaviors 

The QLD shell-structures in Eq. (1561) affect physical values in the super-anisotropic sys- 
tems. Here we show the Fermi energy ep for an example. The Fermi energy is determined 
by the particle number condition, 

N =7^ f dr f d Pf( r 'P) (61) 



(2vr) 3 d 

with the Wigner function in Eq. (|60|) . In particular, we obtain 



W-'d!(n^Ar) ' (62) 
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in the d- dimensional TEA with the zero-point energy shift, = lu c , e = u> z /2, and 
e[, 3) = 0, according to Eqs. (p]), (1551). and (j55) (generally Eq. (jB7j) ). 

In Fig. [3X we show the Fermi energies in the oblate and prolate deformed potentials when 
N = 10 4 and 10 3 , respectively. When the anisotropy of the trap is not so large, the exact 
values of the Fermi energy (the open circles) approximately agrees with that in 3DTFA (the 
short-dashed lines); as the anisotropy increases, it shows crossover behavior from that in 
3DTFA to 2DTFA or 1DTFA and exactly agrees with that in 2DTFA (the long-dashed line 
in Fig. [3^) when ep = (3/2)co z and N = {uj z /uj c ) 2 /2 in the oblate deformed potentials and 
in 1DTFA (the long-dashed line in Fig. [3b) when = 2u c and N = u c / u z in the prolate 
deformed potentials . Here one should note the results in ETFA (the solid lines), which 
include the QLD shell-structures only below e = (3/2)u) z or e c = 2uj c , i.e. the lowest shell 
state (q z = or q c = 0), as shown in Eq. (1601) . These results roughly reproduce the global 
crossover behaviors and slightly different from the exact calculations in the middle range 
of the anisotropy owing to the effect of the upper QLD shell-structures with the hierarchic 
structures. 

As shown in the behavior of the Fermi energy in Fig. [31 the QLD shell-structures induce 
the crossover behavior between the 3DTF and 2/1DTF gases. The crossover behavior is 
an important feature of the trapped quantum gases. The phase diagrams of the weekly 
interacting Fermi gases can be determined by the ratio of ep and u) Z)C (^> to c ,z) according to 
the dimensional analysis. We briefly describe it in Fig. [H 

B. Determination of the parameters 

Using ETFA, we now decide the parameters needed in the scaling and sum-rule methods. 
According to the perturbation theory, we expand these parameters by the coupling constant 
g and take the zero-th and first order terms: e.g. Kj ~ Kj + Kj\ and Uj ~ Ujo + Uji for 
j = c and z. 

According to the Wigner functions in Eq. f[56|) . the zero-th order terms of the potential 
energies in the c and z-direction motion, U c0 and U z0 , become 




(63) 
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FIG. 3: (color online) The Fermi energies in the oblate deformed potentials when N = 10 4 (a) and 
in the prolate deformed potentials when iV = 10 3 (b). The open circles represent the exact results; 
the solid, short-dashed, and long-dashed lines represent the results in ETFA ( e c = (3/2)u; 2 (a) or 
2uj c (b)), 3DTFA, and 2DTFA (a) or 1DTFA (b), respectively. 

^o=E 6 J/ Q(e F -e 2 ) (64) 

n z =0 c 

in the oblate deformed potentials and 

OO CO / \ 

n x =0 n y =0 " 
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FIG. 4: (color online) The phase diagrams of the weekly interacting Fermi gases in the anisotropic 
traps. 
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\2 



n x =0n y =0 

in the prolate deformed potentials, where we use the e z and e c = e x + e y in Eq. (j4~3l) and the 
integral formulas of the Hermite functions. By smoothing the residual sums in Eqs. (I63l) - (l66|h 
Y^m ~^ I dej/uj (j = x,y,z), we obtain the results in 3DTFA: 

U<$ = 2US = \{^ z ulN^I\ (67) 

By taking n z = and n x = n y = in the residual sums in Eqs. fl63l) - fl66l) . we obtain the 
results in 2DTFA and 1DTFA, respectively: 

U$ = ^NV\ U$ = ^N (68) 

and 

o2> = >, U% = ^N>. (69) 
In Fig. we show the zero-th order terms of the potential energies U c o and U z q in 

Eqs. f l63|) - fl66|) in the oblate deformed potentials when N = 10 4 and in the prolate deformed 

potentials when N = 10 3 , respectively. The crossover behavior and the QLD shell-structures 

also appear in these plots. 

The mass parameters in Eq. ( 13TT) is determined by the potential energies owing to the 

definitions: 

M c , z = -^rU c>z . (70) 



i > 2 
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FIG. 5: (color online) The potential energies, U c o [(a) (c)] and U z q [(b) (d)], in the oblate deformed 
potentials when N = 10 4 [(a) (b)] and in the prolate deformed potentials when N = 10 3 [(c) 
(d)]. The open circles represent the exact results; the solid, short-dashed, and long-dashed lines 
represent the results in ETFA (e c = (3/2) uj z [(a) (b)] or 2uj c [(c) (d)]), 3DTFA, and 2DTFA [(a) 
(b)] or 1DTFA [(c) (d)], respectively. 



The zero-th order terms of the kinetic energies equal those of the potential energies 
according to the virial theorem in Eqs. (|33|) and ([Mil : 



for j = c and z. 

In the first order terms, the perturbation theory gives 



(71) 



K n = -U. 



(72) 



for j = c and z owing to the variable separation of the noninteracting ground state; further- 
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more, the virial theorem in Eqs. (1331) and (|34|) gives 

E int = U cl - K cl = 2(U zl - K zl ). (73) 

According to Eqs. (1721) and (I731 . we obtain the first order terms: 

U cl = -K cl = ^E int , U zl = -K zl = -E iQt . (74) 

The interaction energy E- in % in the crossover process varies between those in the d- 
dimensional TFA, 

E$=iptogyfe%Ni, (75) 

where (3^ = 128^/(945^) w 0.00756, 0® = 1/(24tt 3 / 2 ) w 0.00748, and /?W = 
V2/(6n 3 ) ta 0.00760. Owing to /? (3) w /5 (2) w the £$2 in Eq. (ESI) are nearly same 
values, so that the interaction energy has little influence on the QLD shell-structures, and 
we approximate it as 

E^t « (76) 

in the following sections. 



IV. BEHAVIORS OF THE BREATHING OSCILLATIONS 

We here show the calculational results by using the parameters determined in the previous 
section. Firstly we calculate the normal mode frequencies of the breathing oscillations in 
the scaling and sum-rule methods. Next we demonstrate the time-evolution in the present 
framework. 



A. Normal mode frequencies 

According to the eigen-equations (Eq. fl36l) ) with Eqs. fl3T|) and (1321 . we obtain the normal 
mode frequencies of the breathing oscillations, fi CjZ , up to the first order of the interaction 
energy as 




n > - w-Tir + a i Mj(Kj + m fl " + ° (EL) (77) 
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for j = c and z, where a c = 1 and a z = 1/4. These solutions are also expanded by the 
coupling constant g: Qj ~ f^o + Qji, where 



(K j0 + U j0) 



^^r-n^- (78) 



and 

Jl " 2 V^'o + ^o M j0 K j0 + U j0 J- 1 J 

The zero-th order terms correspond to those of the noninteracting gases. By substituting 

Eqs. (COD and ([71]) into Eq. (JTSj), we obtain 

fi i0 = 2to>j. (80) 

This result only depends on the trap frequencies and does not be affected by the QLD 
shell-structures. 

The first order terms are also obtained as 

nci = o, n, 1 = -^^ (si) 

by substituting Eqs. flTOj) . (J7T|) . and f[71"|) into Eq. f[79"j) . This result reveals that the interaction 
and QLD shell-structures affect only the frequency of the longitudinal normal mode, but not 
that of the transverse normal mode. The vanishment in the transverse normal mode reflects 
a typical phenomenon of the transverse mode because the mixing, i.e. the off-diagonal 
elements in Eq. (1321) . has no influence on the first order terms in Eq. (I8ip . 

According to Eq. ( 1HT1) . the QLD shell-structures contribute the frequency of the longitu- 
dinal normal mode principally through the potential energy U z o because Eq. (1811) contains 
only the two parameters, U z q and E- m t, and the QLD shell-structures have little influence 
on the interaction energy E int as shown in Eq. (176]) . Thus the frequency of the longitudinal 
normal mode exhibits the crossover behaviors associated with those of the potential energy 
U z q as shown in Fig. [5j At the (^-dimensional TEA limits (Eqs. (j67j) - (j69"l) ). the corresponding 
frequency becomes 

( 2) = _lEto 

zl 2 N K ' 



and 



nff = -kw- (84) 
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In Fig. [6j we show the first order term in the frequency of the longitudinal normal mode, 
Q z i, in Eq. ([51]) by using the parameters in Eqs. flM|) . fl66|) . and fl76|) . Here one can see the 
crossover behaviors in the exact results (the open circles) between the results in 3DTFA (the 
short-dashed lines) and 2DTFA (a) or 1DTFA (b) (the long-dashed lines). 
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FIG. 6: (color online) The first order term in the longitudinal normal mode frequency, Q z i, in 
the oblate deformed potentials when N = 10 4 (a) and in the prolate deformed potentials when 
N = 10 3 (b). The open circles represent the exact results. The short-dashed and long-dashed lines 
represent the results in 3DTFA and 2DTFA (a) or 1DTFA (b), respectively. 

In Fig. we show comparison between the exact results in Fig. |6] (the open circles) and 
the results in ETFA at e c = (3/2)u z , (5/2)a> 2 , and (7/2)ou z (a) or 2u c , 3u c , and 4u c (b) 
(the short-dashed, long-dashed, and solid lines, respectively). There appear three ranges 
reflecting the hierarchic structures: (i) uj z /oj c > \/2N (a) or uj c /uj z > N (b). (ii) \/2N > 
uo z /uj c > y/(2/5)N (a) or > co c /uj z > N/4 (b). (iii) ^{2/h)N > u z /u c ( a) or N/4 > u) c /uj z 
(b). 

In range (i), the exact results agree with those in 2DTFA (a) or 1DTFA (b), so that all 
results in ETFA reproduce the exact values. 

In range (ii), the QLD shell-structures below the first shell-structure only contribute 
to the exact results. Thus the results in ETFA at e c > (5/2)u z (a) or 3u c (b) can only 
reproduce the exact values. The difference between the results in exact calculations and 
ETFA at e c = (3/2)u z (a) or 2u c (b) indicates the effect of the QLD shell-structures of the 
first shell-structure. 

In range (iii), the QLD shell-structures below the second shell-structure only contribute 
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FIG. 7: (color online) Same as Fig. [H but the short-dashed, long-dashed, and solid lines represent 
the results in ETFA at e c = (3/2)u z , (5/2)u z , and (7/2)uj z (a) or 2u c , 3uj c , and 4w c (b), respectively. 

to the results. Thus the results in ETFA at e c > (7/2)u z (a) or 5u c (b) can only reproduce 
the exact values. There are additional detailed structures associated with the hierarchic 
structures; however it is difficult to find them up in Fig. because the difference between 
the results in exact calculations and ETFA is visually small. 

Finally we comment on the particle number dependence on the above results. When the 
particle number varies, there principally appear two kinds of quantitative change. One is the 
change of the values at the 2DTF and 1DTF limit as shown in Eqs. (IS2l -(1541): ftg 5 oc N 1/6 , 
Q z 2 i oc N 1 / 2 , and fir^i oc N^ 1 ^ 2 . The other is the change of the critical anisotropy points 
between the ranges (i) and (ii), which are proportional to N 1 ^ 2 for the oblate deformed 
potentials and to N for the prolate deformed potentials. 

B. Time-evolution 



We here show the results of the time-evolution in the present framework by using the 
parameters determined in the previous section. 

In the scaling method, the collective coordinates X c ,z{t) determine the behaviors of the 
breathing oscillations as described in Eq. (|T6|) . In order to observe the collective coordinates, 
it is useful to measure the projected mean square radius s C)Z (t) defined in Eq. (|T3l) . which 
are proportional to X c>z {t). We here demonstrate the time-evolution of s CjZ (t) (oc A,. - (I)) 
by solving the dynamical system (Eqs. ( |22l) and ( |30l) ) in the classical mechanics 36(; the 
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equations of motion, up to the first order of the the interaction energy E- mt , become 

d ^ = -^AK-icK) (85) 

and 



dt 
with 



f/ "'' V " - -n 2 z (X z - lz X c ) (86) 



E\ r 



lc, Z = ~Q2~J]~- ( 87 ) 

In Fig. [HI we show the time evolution of the collective coordinates of the transverse and 
longitudinal oscillations, s c (t) and s z (t), in the oblate deformed potential (u z /u c = 80) 
when N = 10 4 and in the prolate deformed potential {uo c /uo z = 600) when N = 10 3 . Here we 
choose an initial condition, s c (t = 0) = s z (t = 0) = so and s c (t = 0) = s z (t = 0) = 0, which 
correspond to an experimental situation by sudden change of the trap frequency at t = 0. 
The short-dashed lines represent the results of the non-interacting gases (g — 0); the solid 
lines and the dotted separations represent the results and the periods of the normal modes 
(Eqs. (|SQD and flUD) when g = 0.3^ 1/2 for the oblate case and g = 0.2u z for the prolate 
case, respectively 

As shown in Fig. [HI the inter-atomic interaction principally induces two kinds of influence: 
frequency decrement in the longitudinal oscillations (Fig. [8b and [HU) and shift of the center 
position of the faster oscillations (Fig. [8b and [Hb). 

The frequency decrements are confirmed from the agreements of the frequencies in the 
results and the corresponding normal mode frequencies as shown in Fig. [8] (the solid lines 
and dotted separations) and Eqs. ( J85|) and ( |86j) : according to Eq. ( IHTi) . the inter-atomic 
interaction keeps the normal mode frequencies in the transverse oscillations, and reduces 
(increases) those in the longitudinal oscillations when g > (< 0), so that the frequency 
decrements occur only in the longitudinal oscillations. 

The center shifts in the faster oscillations are induced by the mixing of the faster and 
slower oscillations. The mixing effects of the faster oscillations on the slower oscillations 
vanish owing to the coarse graining in the long time scales corresponding to the periods of 
the slower oscillations 37| ; however the mixing effects of the slower oscillations on the faster 



oscillations remain as the center shifts, which are nearly constant in the small time scales 
and exhibit beat structures in the long time scales. The center shifts can also be seen in 
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FIG. 8: (color online) The time evolution of the collective coordinates s c [(a) (c)] and s z [(b) (d)] 
in the oblate deformed potential (uj z /u> c = 80) when iV = 10 4 [(a) (b)] and in the prolate deformed 
potential (uj c /oj z = 600) when N = 10 3 [(c) (d)]. The short-dashed lines represent the results when 
5 = 0; the solid lines and the dotted separations represent the results and the periods of the normal 
modes when g = 0.3u> c [(a) (b)] or 0.2u z [(c) (d)], respectively. 

Eqs. ( 1851) and (1561) : when the X z in Eq. (I85I) (A c in Eq. (!86|) ) is nearly constant in time as 
the slower oscillations, the mixing term 7 C A^ in Eq. (185]) (7^A C in Eq. (186]) ) gives the center 
shift. 

Here it should be noted that the mixing does not affect the frequency of the normal modes 
in the first order terms, so that the mixing effects principally appear in the center shifts of 
the faster oscillations; furthermore we can say that the faster oscillations can exhibit well 
collectivity at the anisotropic limits because the atoms there take same one-particle energy- 
level in the narrower direction motions, and the one-particle wave-functions only within the 
energy- level principally determine the behaviors of the faster oscillations just like as the 
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collective oscillations of the BECs, where the strength functions take nearly a single peak in 
the excited energy. Thus we point out that the faster oscillations must include the detailed 
informations of the super-anisotropic deformed gases. 



V. SUMMARY AND OUTLOOK 

In the present work, we study the breathing oscillations of the fermion condensates in 
the anisotropic harmonic oscillator traps. If anisotropy of the traps is not so large, the 
gases behave as the 3D gases; otherwise the gases exhibit the QLD properties owing to the 
QLD shell-structures associated with the trap anisotropy. We theoretically reveal the effects 
of the QLD shell-structures on the breathing oscillations of the week-interacting fermion 
condensates with same number densities and describe the crossover behaviors between the 
3D and QLD gases by developing TEA to the super-anisotropic systems. 

The collective frequencies of the breathing oscillations are calculated in the scaling and 
sum-rule methods with the parameters determined in the perturbation theory and ETFA up 
to the first order of the inter-atomic interaction. As a result, it turned out that the effects 
of the QLD shell-structures disappear in the frequency of the transverse normal mode and 
appear only in that of the longitudinal normal mode with the hierarchic structures associated 
with the trap anisotropy. The results in ETFA with the finite cut-off energy theoretically 
reproduce the hierarchic structures below the cut-off energy, and agree with the exact results 
for the proper cut-off energy. 

Finally we also demonstrate the time-evolution of the breathing oscillations in the present 
framework. In conclusion, we clarify the two kinds of influence of the inter-atomic interaction 
on the transverse and longitudinal oscillations: the frequency decrements (increments) in 
the longitudinal oscillations when g > (< 0) and the center shifts in the faster oscillations. 

The results in the present work give a theoretical prediction for the precise observations in 
the weekly interacting gases. If the inter-atomic interaction is very strong, our method is not 
applicable, and the oscillation properties may be different from our results. In addition, when 
the oscillations are including more than one strong oscillation modes, e.g. near the phase- 
separation or collapse, the scaling and sum-rule methods in this work should not be applied. 
In those cases, one should do more detailed investigation in the full-microscopic RPA or other 



direct time-evolution approaches compatible with the super-anisotropic systems 32j, [33]. 
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Asymmetric two-component gases in the super-anisotropic traps are another interesting 
topic. In that case, the breathing oscillations implicitly include not only the in-phase oscil- 
lations but also the out-of-phase oscillations, and exhibit coupled oscillations of the in-phase 
and out-of-phase oscillations |l( ■ We should study the ground state properties and collective 
excitations in another paper. 
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APPENDIX A: SCALING AND SUM-RULE METHODS 

The scaling method gives the density matrix in the collective oscillation by an excitation 
operator O: 

g(X,X) = U(X,X)g U\X,X) (Al) 

with a unitary operator 

U(X,X) = e ao e x ^°\ (A2) 

where A is the scaling amplitude of the collective oscillation and A indicates the time deriva- 
tive. Here the density matrices are normalized, tr (g) = tr (g ) = 1, and the expectation 
values of an arbitrary operator A are given by (A) = tr (gA) and \A} Q — tr (go A). 
The time- dependent variational principle, 

5 J dt(ij t - H) =S f dt £[A,A]=0, (A3) 

gives the Euler-Lagrange equation for the generalized dynamical variable A. By substituting 
Eq. flAlj) into Eq. flA3j) . we can write the Lagrangian C in Eq. flA3j) as 

C[X, A] ~ ([[O, H] , 0]) O A 2 - l([[0, H] , [H, [H, 0}}]\X* (A4) 

to the second order in A beside the constant term of A. 

When the density matrices are pure states, i.e. £o = | x I / o)( x I / o|, the coefficients in Eq. (1A4I) 
becomes 

([[O 1 ff\,O)) = M 1 , ±{[[O,H\,[H,[H,O]))) = Ms (A5) 
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with the n-th moment 

M n = J2(E U - E ) n |(*„|0|tf o )| 2 , (A6) 

where \ is the u-th state with the energy E v in a complete set. As a result, the frequency 
of the collective oscillation can be written as 



This result is equivalent to that in the sum-rule method. 

APPENDIX B: THOMAS-FERMI APPROXIMATION 

In the one-particle picture, the Wigner function in Eq. ( )45j) can be described as 

f(r,p) = Jde F(e)D(s)f^\r,p;s), (Bl) 

where we use the Fermi distribution function F(e) = l/(exp(/3(e — //p)) + 1) and introduce 
the one-particle Wigner function 

/<«*> (r , p; e) = £ 5 £ , en f ds e'^n (r - |) n (r + |) (B2) 
and energy-level density 

D{e) = ^5(e-e n ) (B3) 

with one-particle energy e„. 

In TFA, the Wigner function in Eq. (1B1I) is evaluated with a semi-classical way: 

D(5)/(° p )(r,p;e) «5(e- ec i(r,p)) (B4) 

or 

/(r,p)^F(e cl (r,p)) (B5) 

with the classical energy e c i(r,p), where we keep the quantum statistics, F(e), and approx- 
imate the quantum motion and energy-level density into these of the classical mechanics. 
This result is also obtained by the semi-classical h expansion [341 ] . 

The classical energy-level density can be obtained by integrating Eq. ( 1B4I) in the phase- 
space owing to the normalization of the one-particle Wigner function: e.g. in the d- 
dimensional harmonic oscillator potentials {d — 1, 2, . . . ), 

eci(r,p) = ^ + ^X>i r ;> ( B6 ) 

3=1 
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the classical energy level densities D^ d \e) become 



£ d-i 



d 



-i 



with the Heaviside's step function (x). 



\j=i 



[1] T. Maruyama and T. Nishimura, Phys. Rev. A 75, 033611 (2007). 

[2] For reviews, see A. S. Parkins and H. D. F. Walls, Phys. Rep. 303, 1 (1998); F. Dalfovo, S. 
Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999); W. Ketterle, 
D. S. Durfee, and D. M. Stamper-Kum, in Bose-Einstein Condensation in Atomic Gases, 
Proceedings of International School of Physics "Enrico Fermi", edited by M. Ingusico, S. 
Stringari and C. Wieman (IOS Press, Amsterdam, 1999); C. J. Pethik and H. Smith, Bose- 
Einstein Condensation in Dilute Gases (Cambridge University Press, Cambridge, 2002). 

[3] E. A. Cornell and C. E. Wieman, Rev. Mod. Phys. 74, 875 (2002); W. Ketterle, Rev. Mod. 
Phys. 74, 1131 (2002). 

[4] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Science 
269, 198 (1995); K. B. Davis, M. O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, 
D. M. Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 

[5] B. DeMarco and D. S. Jin, Science 285, 1703 (1999). 

[6] S. R. Granade, M. E. Gehm, K. M. O'Hara, and J. E. Thomas, Phys. Rev. Lett. 88, 120405 
(2002). 

[7] A. G. Tuscott, K. E. Strecker, W. I. McAlexander, G. B. Parridge, and R. G. Hullet, Science 
291, 2570 (2001). 

[8] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, T. Bourdel, J. Cubizolles, and C. Sa- 
lomon, Phys. Rev. Lett. 87, 080403 (2001). 
[9] Z. Hadzibabic, C. A. Stan, K. Dieckmann, S. Gupta, M. W. Zwierlein, A. Gorlitz, and W. 
Ketterle, Phys. Rev. Lett. 88, 160401 (2002). 
[10] Z. Hadzibabic, S. Gupta, C. A. Stan, C. H. Schunck, M. W. Zwierlein, K. Dieckmann, and 

W. Ketterle, Phys. Rev. Lett. 91, 160401 (2003). 
[11] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 92, 040403 (2004); M. Bartenstein, 
A. Altmeyer, S. Riedl, et.al., Phys. Rev. Lett. 92, 120401 (2004); M. N. Zwierlein, C. A. Stan, 

30 



C. H. Schunck, et.al., Phys. Rev. Lett. 92, 120403 (2004). 
[12] For reviews, see S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 80, 1215 
(2008). 

[13] T. Sogo and H. Yabu, Phys. Rev. A 66, 043611 (2002). 

[14] A. Gorlitz, J. M. Vogels, A. E. Leanhardt et.al., Phys. Rev. Lett. 87, 130402 (2001). 

[15] S. Dettmer, D. Hellweg, P. Ryytty et.al., Phys. Rev. Lett. 87, 160406 (2001). 

[16] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, T. Bourdel, J. Cubizolles, and C. Sa- 
lomon, Phys. Rev. Lett. 87, 080403 (2001). 

[17] T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. Esslinger, Phys. Rev. Lett. 92, 130403 
(2004). 

[18] H. Ott, J. Fortagh, G. Schlotterbeck, A. Grossmann, and C. Zimmermann, Phys. Rev. Lett. 
87, 230401 (2001). 

[19] W. Hansel, P. Hommelhoff, T. W. Hansch, J. Reichel, Nature 413, 498 (2001). 

[20] J. H. Thywissen, M. Olshanii, G. Zabow et.al., Eur. Phys. J. D 7, 361 (1999). 

[21] S. Schneider, A. Kasper, Ch. vom Hagen et.al., Phys. Rev. A 67, 023612 (2003). 

[22] For review, see L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Clarendon Press, 

Oxford, 2003); D. S. Petrov, D. M. Gangardt, and G. V. Shlyapnikov, J. Phys. IV 116, 3 

(2006). 

[23] G. M. Bruun, Phys. Rev. A 63, 043408 (2001). 

[24] K. Goral, M. Brewczyk, and K. Rzazewski, Phys. Rev. A 67, 025601 (2003). 
[25] G. M. Bruun and B. R. Mottelson, Phys. Rev. Lett 87, 270403 (2001). 
[26] L. Vichi and S. Stringari, Phys. Rev. A 60, 4734 (1999). 

[27] G. F. Bertsch, Nucl. Phys. A 249, 253 (1975); G. F. Bertsch and K. Strieker, Phys. Rev. C 
13, 1312 (1976); D. M. Brink and R. Leobardi, Nucl. Phys. A 258, 285 (1976); T. Suzuki, 
Prog. Theor. Phys. 64, 1627 (1980). 

[28] O. Bohigas, A. M. Lane, and J. Martorell, Phys. Rep. 51, 267 (1979). 

[29] T. Maruyama and G. F. Bertsch, Phys. Rev. A 73, 013610 (2006). 

[30] A. L. Fetter and J. D. Walecka, Quantum theory of many-particle systems (McGrawHill, NY, 
1971). 

[31] The precise structures of the one-particle energy-level densities are well-studied in the Nilsson 
model for nuclear physics; See P. Ring and P. Schuck, The Nuclear Many-Body Problem 

31 



(Springer, New York, 2000). 
[32] T. Maruyama, H. Yabu, and T. Suzuki, Phys. Rev. A 72, 013609 (2005). 
[33] T. Maruyama and G. F. Bertsch, Phys. Rev. A 77, 063611 (2008). 

[34] R. G. Parr and W. Yang Weitao, Density-functional Theory of Atoms And Molecules (Oxford 
University Press, New York, 1989). 

[35] If we consider the quantum number of the coarse mode as a component number, the 3D gas 
can be considered as a multi-component mixture of the QLD gases. Thus the dimensionality 
is related to the multiplicity. 

[36] The X c ,z(t) are taken as the c- numbers in the scaling method, so that the classical system 
in Eqs. (|22|) and (|30p must contain the quantum fluctuations as mean-values. This property 
originally comes from the time-dependent HF theory. 

[37] The coarse graining is due to the multi-time-scale structures derived from the QLD shell- 
structures, and yields the QLD properties in the slowest time-scale with renormalization of 
the freedom degrees associated with the other time-scales as plotted in Fig. 



32 



